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PREDICTED  VALUES  OF  THE  VISCOSITY  AND  THERMAL  CONDUCTIVITY 
COEFFICIENTS  OF  NITROUS  OXIDE" 

Howard  J.  M.  Hanley 

The  viscosity  and  thermal  conductivity  coefficients  of  nitrous 
oxide  are  calculated  for  temperatures  between  180  and  900  K  (330  to 
1600°R)  for  pressures  to  23  MPa  (^  3500  psi)  .   Tables  of  values  are 
presented.   Two  mixtures  with  carbon  dioxide  are  also  discussed.   These 
transport  coefficients  (for  the  pure  fluid  and  for  the  mixtures)  were 
predicted  from  thermodynamic  data.   Details  of  the  prediction  procedure 
are  presented.   Estimates  of  the  accuracy  of  the  tabular  values  are 
+  6%  for  the  viscosity  and  +  8%  for  the  thermal  conductivity. 

Key  words:   Carbon  dioxide;  corresponding  states;  mixtures;  nitrous  oxide; 
prediction;  thermal  conductivity;  transport  property;  viscosity. 

1.   INTRODUCTION 

In  this  report  we  discuss  the  transport  coefficients  —  the  viscosity 
coefficient  (f))  and  the  thermal  conductivity  coefficient  (A)  —  of  nitrous  oxide, 
N~0.   Tables  of  values  are  presented  from  0.1  to  23  MPa  (approximately  15  to 
3500  psi)  for  temperatures  between  180  and  900  K  (approximately  330  to  1600°R) . 
Two  selected' mixtures  of  nitrous  oxide  and  carbon  dioxide  are  considered  briefly, 

Nitrous  oxide  is  a  common  substance  yet  the  transport  coefficients  have  not 

been  measured  over  a  wide  temperature  and  pressure  range.   In  fact,  the  data 

are  so  scarce  we  prefer  to  obtain  the  coefficients  by  calculation  [1]:  the 
limited  data  are  used  only  to  check  the  procedure  where  possible. 

2.   DATA  SURVEY 

A  bibliography  of  data  sources  for  the  thermodynamic  and  transport 
properties  of  nitrous  oxide  was  prepared  by  the  Cryogenic  Data  Center,  National 
Bureau  of  Standards,  Boulder,  Colorado.   The  appropriate  references  for  the 
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viscosity  and  thermal  conductivity  are  as  follows:   viscosity  coefficient, 
references  [2]-[8];  thermal  conductivity  coefficient,  references  [9]-[23]. 

We  have  correlated  the  transport  coefficients  of  several  simple  fluids 
[Ar,  0„,  N„ ,  CH, ,  C~H,]  in  a  series  of  papers  [24].   Criteria  were  set  up  to 
evaluate  data  critically,  and  an  equation  to  represent  the  coefficients  was 
proposed.   Unfortunately,  this  procedure  cannot  be  applied  to  nitrous  oxide  at 
this  time  because  suitable  data  are  not  available  for  analysis.   We  have  been 
able  to  find  only  two  data  sets  for  the  dense  gas  and  liquid,  those  from 
references  [19]  and  [23],  for  the  thermal  conductivity.   The  other  references 
cover  the  dilute  gas  only. 

As  remarked,  therefore,  we  decided  to  present  tables  based  on  a  prediction 
procedure,  which  will  be  outlined  in  the  next  section. 

3.   PREDICTION  METHOD 

The  method  used  is  an  application  of  a  procedure  introduced  in  detail  in 
reference  [1].   An  outline  is  given  here. 

As  is  well  known,  if  classical  two-parameter  corresponding  states  applies 
to  two  fluids  designated  a  and  o,  one  can  write  the  viscosity  of  one  at  a  density 
(p)  and  temperature  (T) ,  in  terms  of  the  other: 

n  (p,t)  =  n  (p   ,  t   )(-*)    -«   (-«)  (i) 


lpo 


Also  for  the  thermal  conductivity, 


1/2/    c\2/3/^cvl/2 


»(pvv|  (7)  (7) 


Vp.t)  =  up„_.  t„  j\~\   \-a  1^1  (2) 


In  the  above  equations,  M  is  the  molecular  weight,  and  p    and  T    are  the 
density  and  temperature  respectively,  for  fluid   o   to  correspond  to  the  density 
and  temperature  for  fluids  a.   The  variables  can  be  expressed  in  terms  of  the 
critical  parameters,  p°  and  TC;  thus,  p^   =  p(pC/p*~)  and  T    =  T(T  /T  ).   The 

CXjO  OCX  \X  y  O         (JUL 


coefficients  r|  and  A   also  have  to  be  scaled  as  shown  in  equations  (1)  and  (2). 
o      o 

It  is  common  practise  to  calculate  the  coefficients  for  a  given  those  for 
o  —  which  is  then  treated  as  a  reference  fluid. 

Equations  (1)  and  (2)  can  be  extended  for  a  mixture,  x,  if  it  is  assumed 
that  the  mixture  can  be  represented  as  an  equivalent  pure  fluid  —  the  one-fluid 
approximation  —  and  that  the  mixture  is  "conformal,"  that  is,  all  components 
obey  the  two-parameter  corresponding  states  law.   Hence  one  can  write 

1/2/  c\2/3/_c\l/2 


,«(?) 


V?>T)  "  Vpa,o-  Ta,o'  \tl     \f>     ; 


and  similarly  for  the  thermal  conductivity,  A  .  Equation  (3)  can  be  derived 
under  well-defined  assumptions  [25] .  The  derivation  further  yields  a  set  of 
mixing  rules: 


Ct  p 


and 


TC  (pC)  1     =      Z  Z  x  x0  TC0  (pC0)  1  (5) 

x  Kx        „  a  MP   ag 

where  x  is  the  mole  fraction  of  a.   The  mixing  rule  used  in  this  work  for  the 
mass  is  M  =  Z  x  M   (which  differs  slightly  from  that  proposed  in  reference  [25]) 
The  a$  variables  of  equations  (4)  and  (5)  are  given  by 

Tc   =  F        (t°   TC  )ly/2  (6) 


and 


/  c  \-l      i 
(paS}    =  ^a3 


1  ,   c  ,-1/3  .  1  .  c  .-1/3 

2  (paa>     +  2  (p33} 


(7) 


where  E,   a      and  \\l   D      are  binary  interaction  parameters,  which  are  best  obtained 

CXp  Ctp 

from  experiment. 

The  mixing  rules,  equations  (6)  and  (7),  are  consistent  with  the  mixing 
rules  derived  for  thermodynamic  properties  using  the  one-fluid  approximation. 
If,  therefore,  the  binary  interaction  parameters  are  assumed  to  be  those  valid 
for  thermodynamic  properties,  one  has  a  procedure  to  predict  the  viscosity  and 
thermal  conductivity  coefficients  of  mixture  x,   given  the  equivalent 
coefficients  of  the  reference  fluid. 

Equations  (l)-(3)  are  not  general  in  that  two-parameter  corresponding 
states  is  invalid  for  polyatomic  fluids,  or  for  mixtures  containing  polyatomic 
molecules.   The  situation  is,  of  course,  paralleled  with  respect  to  the 
thermodynamic  properties  and  has  been  discussed  extensively  in  this  latter 
context.   For  example,  a  recent  approach  due  to  Leland  [26]  and  to  Rowlinson  [27], 
is  of  interest  here.   According  to  Leland  and  Rowlinson  a  third  parameter, 
co  —  which  can  be  taken  as  the  Pitzer  acentric  factor  —  is  introduced,  but  the 
framework  of  simple  two-parameter  corresponding  states  is  preserved,  if  one 
considers  the  scaling  functions 

f      =   (TC  /TC)  0      ;      h         (PC/P°  )  <J>  (8) 

aa,o      aa  o   aa,o       aa,o      o  aa  Taa,o 


for  fluid  a  with  respect  to  o.   The  terms  6     and  6     are  called  shape  factors 

r  aa,o     aa,o 

and  are  weakly  varying  functions  of  temperature  and  density: 

6      =  1  +  (co  -co  )F(T,p);    <J)      =  1  +  (co  -co  )G(T,p)  (9) 

aa,o  aa  o    '   '    Yaa,o  aa  o 

where  co   and  co  are  acentric  factors  for  fluids  a  and  o,  respectively.   Using 

aa      o 

the  shape  factors,  the  compressibility  factor  (for  example)  of  a  would  be  given 

by  the  relation  Z  (p,T)  =  Z  (ph     ,  T/f     ).   With  the  appropriate  expression 
3  a        o   aa,o     aa,o 

for  the  Helmholtz  free  energy,  the  thermodynamic  properties  of  a  are  completely 
defined  in  terms  of  the  properties  of  o. 

It  is  a  logical  step  to  introduce  shape  factors  into  the  viscosity  and 
thermal  conductivity  expressions:   for  fluid  a,  therefore 


n  (p,t)    =    n    (p'     ,  t'    )  FHn  (io) 

a  o       a,o       a,o         a,o 


where 


p'  -      Ph  ;      T'  =      T/f  ;      FHn  =       -£ )         h"2/3   fl/2  (      , 

a,o  aa,o         a,o  aa,o  a,o  \M   /  aa,o     aa,o 


fe)'" 


Similarly, 

A    (p,T)      =      X      (p'       ,    TT      )    FHA  (12) 

a  K'  o       a,o       a,o         a,o 


where 

V 

a,o  \M    /        "aa,o     aa,o 


\M  /  aa,o 


FH*  -     U  h"2/3   £l/2_  (13) 


And  for  mixtures, 

\(P'T)   =  %  (0  Tx,o}  <o  (14) 


A  (p,T)   =   A   (p'   ,  T'   )  FHX  (15) 

x  o   x,o   x,o    x,o 


For  mixing  rules  for  h    and  f    we  take 

x,o  x,o 

h  =ZExx0hD  (16) 

x,o  a  a  3     a3,o 


f  h  =EZxxDf0        hD  (17) 

x,o     x,o  _       a  3     a3,o     a3,o 


with 


1/2 
fa(3,o  "  ^ag  (faa,o  f3g,o)  (18) 


and 


a(3,o     rag 


3l 

;,oJ 


1  ,1/3    1  ,1/3  I3 

2  haa,o  +  2  h*«  -1  (19) 


An  evaluation  of  equations  (10),  (12),  (14)  and  (15)  is  discussed  in 
reference  [1].   In  particular  we  were  interested  if  these  equations  could  be 
used  to  predict  the  transport  coefficients  of  a  or  x,  given  values  for  the 
reference  fluid.   Specifically,  one  might  hope  that  equations  could  be  used 
successfully  if  the  shape  factors,  9  and  (J),  and  the  mixing  interaction 
parameters,  £  R  and  \\)   R,  were  taken  from  a  fit  of  thermodynamic  (PVT)  data. 
Unfortunately,  the  equations  could  not  represent  data  satisfactorily  if  the 
above  constraint  was  imposed. 

We  have,  however,  been  able  to  propose  a  possible  modification;  we  suggested 
that  one  should  consider  expressions  of  the  form 


n  (p,T)   =  n   (p'   ,  T»   )  FHn   Xn    (p,T)  (20) 

a  K'       o   a,o   a,o    a,o  a,o 


and 

A  (p,T)   =   A   (p1   ,  T'   )  FHA   XX    (p,T)  (21) 

a  o   a,o   a,o    a,o  a,o 

with  corresponding  expression  for  the  mixture.   Note  that  equations  (20)  and  (21) 

differ  from  those  introduced  earlier  by  the  factor  X    or  X„   .   Based  on  how 

J  a,o     a,o 

equations  (20)  and  (21)  represent  experiment,  it  turns  out  that  X  should  be 
unity  if  fluids  a  and  o  follow  classical  two-parameter  corresponding  states, 
as  in  the  case  of  the  rare  gases  and  their  mixtures  [25],  but  should  be  a 
function  of  temperature  and  density  otherwise. 

No  formal  equation  is  available  for  X  but  a  possible  expression  was 
determined  in  reference  [1] .   The  expressions  are  based  on  a  study  of  the 
Modified  Enskog  Theory.   This  theory  has  been  discussed  in  detail  in  our 


previous  work  [28].   Following  the  arguments  that  are  given  in  full  in 
reference  [1],  we  obtained  for  the  pure  fluid  (with  corresponding  expressions 
for  the  mixture) 

X^   (p,T)   =   Q    G^     :   X*   (p,T)   =   Q    G*  (22) 

ot,o  ot,o   ot,o      a,o  0L,O      0L,0 


Q    is  function  of  p  and  T  defined  as 
a,o 

0     =   [(bp)  /(bp)n][1-1/exP  (pC/p)3]  (23) 

where  b   is  a  term  given  by  b  =  B  +  TdB/dT,  with  B   the  equilibrium  second 

virial  coefficient.   G    and  G    are  ratios  given,  respectively,  by 

-,a,o,     a,o  b 

Gn   =  []n/[]n  and  []/[].   The  bracket  expressions  are  [28] 
a,o     a   o      a   o 

[]n  =   [1/bpx  +  0.8  +  0.761  bpx] 

(24) 
[]A  =   [1/bpx  +  1.2  +  0.755  bpX] 


with 


PR  \3TJp 


!;^  -  i  <25> 


for  fluids  a  or  o.   R  is  the  gas  constant. 

Notice  that  the  calculation  of  the  correction  factors  X  and  X  require  only 

thermodynamic  information  for  fluid  a  or  x  with  respect  to  the  reference  fluid; 

for  example,  b   [Equation  (23)]  can  be  obtained  from  the  second  virial 

coefficient  of  the  reference  fluid,  B  ,  according  to  the  expression 

o 

(h     B  ) 

b  (T)   =   h     B  +  T  d    aa'°   -  (26) 

a        aa,o  o  dT 


and  similarly  for  the  term  bpx  [equation  (25)]. 


3 . 1   Summary 

To  summarize  the  calculation  procedure.   To  calculate  the  coefficients  of 
a  fluid,  a,  or  mixture,  x,  one  needs  (a)  the  viscosity  and  thermal  conductivity 
coefficients  of  a  reference  fluid,  (b)  the  shape  factors  and,  if  necessary,  the 
mixing  parameters  of  the  fluid  with  respect  to  the  reference  fluid;  these 
quantities  to  be  obtained  from  a  fit  of  thermodynamic  data.   The  viscosity  and 
thermal  conductivity  of  a  or  x  follow  from  the  equations 


n  (P,T)   =  n   (p1   ,  T'   )  FHn   xn  (20) 

aVK'       'o  Ka,o   a,o    a,o  a,o  v  ' 


\    (p,T)   =   A   (p    ,  T?   )  FHA    XA  (21) 

a  o   a,o   a,o    a,o  a,o 


n  (p,t)  =  n  (p'  ,  t'  )  FHn   xn  (27) 

x  o   x,o   x,o    x,o   x,o 


X  (p,T)   =  X      (p'   ,  T'   )  FHA   XA  (28) 

x  o   x,o   x,o    x,o   x,o 


In  this  work  a  is,  of  course,  nitrous  oxide.   Methane  was  chosen  as  the 
reference  fluid  and  numerical  values  of  the  transport  coefficients  were  those 
of  the  correlation  of  Hanley,  Haynes  and  McCarty  [24].   The  equation  of  state 
for  methane  used  here  is  due  to  Goodwin  [29]. 

Shape  factors  for  nitrous  oxide  with  respect  to  methane,  equation  (8), 
were  obtained  by  analyzing  the  PVT  data  of  Couch,  Kobe  and  Hirth  [30].  The 
procedure  is  due  to  McCarty  [31]  and  will  not  be  described  here. 

3.2   Kinetic  Theory 

For  consistency,  the  transport  coefficients  of  nitrous  oxide  for  all 
pressures  and  temperatures  will  be  calculated  by  the  method  summarized  in 
section  3.1.   However,  alternative  approaChs  are  possible  for  the  dilute  gas 
since  some  data  are  available.   In  previous  work  on  the  dilute  gas  [32-34],  we 


have  used  an  empirical  fit  of  data  and/or  calculated  values  via  statistical 
mechanics  and  kinetic  theory.   This  latter  approach  has  several  advantages, 
one  of  which  is  that  statistical  mechanics  and  kinetic  theory  allow  an  assess- 
ment to  be  made  of  systematic  error  by  comparing  independently  measured 
properties  through  the  pair  potential  function.   For  example,  in  references  [32-34], 
we  showed  the  transport  properties  of  the  rare  gases,  and  oxygen  and  nitrogen 
were  consistent  with  their  respective  second  virial  coefficients.   As  a  matter 
of  interest,  nitrous  oxide  is  briefly  discussed  along  these  lines  in  the  Appendix. 

3.3   The  Behavior  of  the  Thermal  Conductivity  Coefficient  in  the  Critical 
Region 

It  is  now  well  known  that  the  thermal  conductivity  behaves  anomalously  in 
the  vicinity  of  the  critical  point:   in  fact  the  coefficient  will  approach 
infinity  as  the  density  and  temperature  approach  their  critical  values.   An 
account  of  this  phenomenon  is  given  in  reference  [35],  and  it  is  considered  in 
our  previous  data  correlations  of  argon,  oxygen,  nitrogen  and  methane  [24]. 
Undoubtedly,  nitrous  oxide  will  also  display  the  anomaly. 

The  anomaly  can  be  calculated  directly  and  is  discussed  in  reference  [35] 
but  it  is  included  indirectly  in  our  corresponding  states  procedure  as  follows. 

The  thermal  conductivity  is  given  by  equation  (21) 

A  (p,T)   =   A   (p    ,  T    )  FHX   XX    (p,T)  (21) 

a  o   a,o   a,o    a,o  a,o 

where  a  E  NO  and   o  E  CH , ,  but  A   can  be  separated  into  two  terms,  the 

conductivity  in  the  absence  of  a  critical  point  anomaly,  A',  and  the  anomaly 

itself,  A  : 
o 

A   =  A'   +  Ac  (29) 

o      o      o 


Hence 


A   =   A'  FHA   XA    +  AC  FHA   XA  (30) 

a      o   a,o  a,o      o   a,o  a,o 


n 

Since  we  have  incorporated  A   in  the  methane  correlation,  A  will  have  an 

r        o  a 

equivalent  contribution. 

4.   COMPARISONS  WITH  DATA 

Comparisons  between  our  calculated  values  of  the  transport  properties  of 
nitrous  oxide  and  data  are  restricted  to  the  dilute  gas  region,  and  to  one  data 
set  for  the  thermal  conductivity  [19].   First  the  dilute  gas. 

We  attempted  to  assess  the  accuracy  of  the  dilute  gas  viscosity  data  using 
the  arguments  discussed  in  reference  [24].   Data  from  references  [6]  and  [7] 
were  selected  as  the  most  reliable  with  an  estimated  accuracy  of  +  3%. 

The  dilute  gas  thermal  conductivity  data  were  also  evaluated  as  far  as 
possible.   Data  from  references  [10],  [14]  and  [18]  were  selected  as  the  more 
reliable  with  an  estimated  accuracy  of  +  6%. 

Following  the  procedure  of  our  previous  work  [24],  the  selected  viscosity 
and  thermal  conductivity  data  were  fitted  to  an  empirical  function  for 
computation  convenience 


n°   =   GV(1)T  X  +   GV(2)T  2/3  +  GV(3)T  1/3  +  GV(4)  +  GV(5)T1/3 
+  GV(6)T2/3  +  GV(7)T  +  GV(8)T4/3  +  GV(9)T5/3 


(31) 


and  similarly  for  A°,  but  with  coefficients  GT(i)(i  =  1...9)  replacing  GV(i) 
in  equation  (31).   n   and  X   refer  to  the  dilute  gas  viscosity  and  thermal 
conductivity  coefficients,  respectively.   Values  of  GV  and  GT  are  listed  in 
table  1. 
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Table  1.   Values  of  the  coefficients  of  equation  (31)  for  the  dilute  gas 
viscosity,  n,  ,  and  the  dilute  gas  thermal  conductivity,  A  .   Units:   temperature 
in  K,  n.   in  ug/(cm*s)  and  A   in  mW/(m*K). 


GV(1 
GV(2 
GV(3 
GV(4 
GV(5 
GV(6 
GV(7 
GV(8 
GV(9 


GT(1 
GT(2 
GT(3 
GT(4 
GT(5 
GT(6 
GT(7 
GT(8 
GT(9 


.2546211337E+07 

-.1670382040E+07 

.2898185505E+06 

.3474760268E+05 

-.1674096936E+05 

.1445438546E+04 

.8767956930E+02 

-.1932828150E+02 

.7959594290E+00 


-.1550251226E+07 

•9682908726E+06 
-.2120883669E+06 

.1779303821E+05 
-.9793550977E+03 

•3480491811E+03 
-.6424000898E+02 

.4582750314E+01 
-.1120206739E+00 
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Figure  1  displays  a  deviation  plot  between  the  experimental  data  fitted  to 
equation  (31),  and  our  estimates  from  equations  (20)  and  (21),  with  the  pressure 
set  at  one  atmosphere.   Our  representation  of  the  dilute  gas  viscosity  coefficient 
can  be  considered  satisfactory,  the  thermal  conductivity  less  so.   for  both 
coefficients  systematic  differences  are  observed  but  it  is  not  clear  whether 
they  are  caused  by  errors  in  the  procedure,  or  are  due  to  errors  in  the  data.   A 
possible  reason  for  the  deviations  in  the  thermal  conductivity  is  that  our 
method  does  not  take  into  account  correctly  the  contribution  of  the  internal 
degrees  of  freedom  of  nitrous  oxide.   This  problem  is  under  further  investigation. 

We  next  compared  our  predicted  conductivities  for  nitrous  oxide  with  the 
dense  gas  and  liquid  data  of  reference  [19].   The  result  is  shown  in  figure  2. 
It  is  seen  that  the  calculated  values  are  within  about  5%,  or  less,  of  the 
measured  values. 

5 .   TABLES 

Tables  of  values  for  the  viscosity  and  thermal  conductivity  coefficients 
of  nitrous  oxide  were  generated  from  equations  (20)  and  (21) .   Tables  2  and  3 
give  the  results  in  the  SI  system  of  units,  tables  4  and  5  in  engineering  units. 
For  convenience,  saturated  liquid  values  are  listed  separately  as  table  6  . 

5.1  Assessment  of  Accuracy 

An  estimation  of  the  accuracy  of  the  tabular  values  is  clearly  difficult 
and  has  to  be  based  mainly  on  how  well  the  prediction  method  can  represent 
other  fluids  for  which  reliable  data  are  available.   For  nitrogen,  ethane, 
propane,  butane  and  carbon  dioxide,  and  their  mixtures,  we  were  generally 
able  to  predict  the  coefficients  to  within  about  6%  using  methane  as  the 
reference  fluid.   Here  we  have  represented  the  dilute  gas  viscosity  data  to  3% 
and  the  dilute  gas  thermal  conductivity  data  to  10%.   The  dense  gas  thermal 


* 

The  number  of  significant  figures  for  an  entry  are  more  than  the  accuracy  of 

the  calculations  warrant  (see  Section  5.1).   Extra  figures  are  included  to 
facilitate  interpolation  if  required. 


12 


conductivity  measurements  of  reference  [19]  have  been  represented  to  within  6%. 
Based  on  the  above  two  factors,  and  our  experience  in  evaluating  data  for  fluids 
similar  to  nitrous  oxide,  we  estimate  the  viscosity  to  be  accurate  to  +  6%  and 
the  thermal  conductivity  to  be  accurate  to  +  8%.   In  the  critical  region  this 
latter  figure  should  be  increased  to  +  20%.   As  in  our  previous  work,  these 
assessments  refer  to  an  accuracy  estimation  on  a  2a  basis. 


13 


Figure  1.   Top  curve:   1  atm  viscosity  percent  deviation  plot  between  values 
predicted  from  equation  (20)  and  the  correlation  equation  (31) .   Percent 
deviation  is  defined  as  n[eq.  (20)]  -  n[eq.  (31) ]*100. 0/n [eq.  (31)]. 
Bottom  curve:   1  atm  thermal  conductivity  deviations  between  equation  (21) 
and  equation  (31) . 
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Figure  2.   Dense  gas  and  liquid  isotherms  for  the  thermal  conductivity 
coefficient  calculated  from  equation  (21)  (curves)  compared  with  the  data  of 
reference  [19]  (open  and  filled  circles). 
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-2 

Table  4.   The  viscosity  coefficient  of  nitrous  oxide.   Pressure  in  10   psi 

[i.e.,  pressure  in  psi  is  obtained  by  multiplying  by  100.0],  temperature  in 
degrees  Rankine,  viscosity  coefficient  in  10   lb/(ft.s)  [i.e.,  to  obtain  X]   in 
lb/(ft.s),  multiply  the  entries  by  10   ]. 
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Table  5-   The  thermal  conductivity  coefficient  of  nitrous  oxide. 

-2 
Pressure  in  10    psi  [i.e.,  to  obtain  the  pressure  in  psi  multiply  by  100.0], 

3 
temperature  in  degrees  Rankine,  thermal  conductivity  in  10   BTU/(ft.  hr.  °R) 

_3 
[i.e.,  to  obtain  A  in  BTU/(ft.  hr.  °R) ,  multiply  the  entries  by  10   ]. 
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Table  6.   Saturated  Liquid  Transport  Coefficients  for  NO 


T 

°F 

■127.03 

•  99.99 

•  50.0 

■  40.0 

■  22.0 

•  4.0 
5.0 

14.0 
23.0 
32.0 
41.0 
50.0 


T 
°R 

332.64 
359.68 
409.66 
419.67 
437.67 
455.67 
464.67 
473.67 
482.67 
491.67 
500.67 
509.67 


n 

lb-mass/f f  s 

.25188E-03 
.19866E-03 
.13464E-03 
.12516E-03 
.10998E-03 
.96711E-04 
•90656E-04 
.84921E-04 
.79464E-04 
.74250E-04 
.68980E-04 
.63901E-04 


BTU/hr«ft«°F 

.11312E+00 
.10210E+00 
.84157E-01 
.80840E-01 
.75034E-01 
•69392E-01 
.66617E-01 
.63863E-01 
.61126E-01 
.58405E-01 
.55536E-01 
.52725E-01 


59.0 
68.0 
77.0 
86.0 
95.9 


518.67 
527.67 
536.67 
545.67 
555.57 


58956E-04 
,54078E-04 
.49158E-04 
.44106E-04 
,  38150E-04 


.50018E-01 
.47495E-01 
.45266E-01 
.43498E-01 
.42357E-01 


30 


6 .   MIXTURES 

Equations  (27)  and  (28)  were  used  to  generate  tables  for  two  mixtures  of 
nitrous  oxide  with  carbon  dioxide:   compositions  C0„  (12%)  -  N„0  (88%)  and 
C0„  (6%)  -  N„0  (0.94%),  respectively.   The  results  are  tables  7-10.   Since 
carbon  dioxide  is  so  similar  to  nitrous  oxide,  our  accuracy  assessments  for 
pure  nitrous  oxide  are  probably  valid  for  mixtures.   No  data  are  available  for 
comparisons. 

7.   CONCLUSION 

We  have  predicted  the  transport  coefficients  of  nitrous  oxide,  and  of 
nitrous  oxide/carbon  dioxide  mixtures  using  a  significant  modification  of  the 
corresponding  states  procedure.   Tabular  values  have  been  listed.   Data  available 
for  comparisons  are  scanty.   More,  and  more  reliable,  data  will  be  required  if 
the  accuracy  of  the  tables  is  to  be  improved. 
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APPENDIX.   Application  of  Statistical  Mechanics 

The  dilute  gas  coefficients,  n   and  X  ,  have  been  correlated  by  the  empirical 
equation  (31).   An  alternative,  however,  is  to  consider  the  kinetic  theory, 
statistical  mechanical,  expressions  such  as 

(TTm  kT)1/2 
°   -  —      — 
n   "  16     2  n(2,2)*  (1A) 

to  9. 


where  m   is  the  molecular  mass,   k   is  Boltzmann's  constant  and   a   is  a  length 

parameter,  characteristic  of  the  intermolecular  pair  potential  function,   (J). 

(2  2 )  * 

The  collision  integral,  0,  ,  is  a  function  of  4>   and  the  temperature.   Clearly 

one  requires  an  expression  for  c})  to  use  equation  (1A)  in  practice. 


The  potential  for  a  nonspherical  molecule  such  as  nitrous  oxide  can  be 

written  as  the  sum  of  a  spherical  contribution,  if>   (i.e.,  4>   depends  on  the 

intermolecular  separation  between  molecules  but  not  on  their  relative  orientations) , 

and  a  nonspherical  contribution,  (J)   ,  which  incorporates  the  appropriate 

ns 

electrostatic  moments  and  molecular  polarizabilities .   Hence 


)   +  0  (2A) 

s     ns 


In  previous  work  we  have  represented  the  spherical  part  by  the  m-6-8  potential 
function  [ 32-34] : 


s 


cj,  /e  =   6±§X  UX   d-  _  P-Yfa-8)  M  d6   /a V  d8 
s        m-6   \r  J  m-6    \r  J  \r  J 


where  d  =  r  .  /a,  with  a  and  r  .   being  defined  by  the  conditions  (J)  (o)  =  0  and 
min  mm  s 

<J>  (r  )  =  -  e,  respectively.   The  parameter   m  represents  the  'hardness'  of  the 
repulsive  term  while  Y  represents  the  'strength'  of  the  inverse  eighth  attraction 
term. 
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It  turns  out  that  the  viscosity  coefficient  is  only  weakly  dependent  on 
the  nonspherical  contribution  so  a  fit  of  viscosity  data  with  equation  (1A) 
will  determine  the  potential  parameters  of  (J)  .   Possible  errors  in  making  this 
simplication  are  discussed  in  reference  [34]. 

Nitrous  Oxide.   Viscosity  data,  [6,7],  of  nitrous  oxide  was  fitted  to  equation 
(1A)  using  the  m-6-8  collision  integrals  from  reference  [36],   The  parameters 
thus  determined  are 

m  =   11  Y  =   3 

a  =   3.65  x  10~10  m         r   =   4.068  x  10_1°m  (4A) 

m 

e/k  =  288.0  K 

A  deviation  curve  comparing  calculated  with  experimental  viscosity  coefficients, 
as  represented  by  the  function  (31),  is  shown  as  figure  3a.   The  deviation  is 
systematic  above  300  K,  but  well  within  the  experimental  uncertainty  of  the 
data. 

The  thermal  conductivity  coefficient  was  also  determined  from  the  corresponding 
kinetic  theory  expression.   The  result  agreeded  with  experiment  to  within  the 
estimated  error  of  6%.   Details  will  not  be  given  in  this  report. 

The  Second  Virial  Coefficient  of  Nitrous  Oxide.   It  is  often  instructive  to 

calculate  the  second  virial  coefficient,   B,   of  a  gas  using  the  potential 

found  suitable  for  the  viscosity.   The  calculation  is  worthwhile  in  itself. 

Further,  however,  the  second  virial  statistical  mechanical  expression  and  the 

second  virial  data  are  independent  of  the  viscosity.   The  calculation  is  thus 

a  valuable  check  on  the  potential,  and  on  the  viscosity  measurements  which  led 

to  the  potential  parameters. 

For  a  polyatomic  molecule,  the  nonspherical  contribution  to  equation  (1), 

(b   ,  has  a  significant  effect  on  the  second  virial  coefficient  and  cannot  be 
ns 

neglected.   In  general  one  has 

(J)    =   (j)  (permanent)  +  (J>(induced)  (5A) 
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where  4> (permanent)  represents  interactions  between  the  permanent  moments  e.g. 

dipole(u)-dipole,  quadrupole(0)-quadrupole,  and  dipole-quadrupole  interactions  — 
while  (J) (induced)  represents  the  contributions  caused  by  the  induction  effects 
of  the  permanent  moments. 

Nitrous  oxide  is  somewhat  unusual  in  that  it  possesses  a  small  dipole 
but  a  large  quadrupole.   In  the  brief  discussion  following,  we  will  consider 
these  moments  but  will  exclude  the  effects  of  (J) (induced)  .   The  potential  for 
nitrous  oxide  is  thus 


)  +  <J>(uu)  +  <K00)  +  <KuO)  (6A) 


where  (j)   is  the  potential  (3A)  with  the  parameters  of  equation  (4A) ;   (J)(uu),  4>(00) 
and  4>(u©)   are  the  dipole-dipole,  quadrupole-quadrupole,  and  dipole-quadrupole 
terms,  respectively. 

Writing  the  equation  of  state  as 


=   1  +  Bp  +  .. .  (7A) 


pNkT 


where  N   is  Avogadro's  number,  the  second  virial  coefficient  is  given  formally 
by 


B   = 


|  if  [exp(-e<f>)-l]  dj^1  d^2  (8A) 


where  3  =  1/kT  and  R1  and  R„  denote  the  position  and  angular  configuration  of 
molecule  1  and  2,  respectively. 

Several  techniques  have  been  proposed  to  integrate  equation  (8A)  of  which 
the  most  common  is  the  expansion  procedure  of  Pople  [37].   One  considers  the 
molecules  to  have  point  moments  and  to  interact  according  to  the  coordinate 
system  shown  in  figure  4. 
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I 


X 


Figure  4. 

Equation  (2A)  for  (j)  is  inserted  into  equation  (8A)  and  the  (J)    term  expanded 
in  powers  of  3-   The  integrations  over  angles  are  then  carried  out  to  yield; 


B*(T*)   =   B(kT/£)/b   =  B*    „  +  B   +B^  +  B^  (9A) 

o      m-6-8      y      0      y0 

3 
Here  b  =  2ttN(J  /3  and  B*  ,  n  is  the  conventional  reduced  second  virial  coefficient 
o  m-6-8 

for  the  spherical  potential,  as  tabulated  in  reference  [36].   For  the  model  shown 
in  figure  4  one  has  for  the  dipole-dipole  potential; 

2 

(J)  =      -  ^r-      [2   Cose,    Cos0o   -   SinG,    Sin6.    Cos(j>]  (10A) 

Tyu  3  12  1  2 

r 


For   (j)00: 


2 

j)Q0    =   |  ~    [i  -  5  cos2e1  -  5  cos2e2  -  15  cos2e1  cos2e2 

r 


(HA) 


+  2  (Sine      Sine      Costj)   -   4   Cos6      Cos62)    ] 


And   for   6   ,-.: 
ry0 


y0    [(Cos60   -   Cose, )  (3   Cos6.    Cos60   -   2    Sin61    Sin6      Cost))) 


y0  „   4      M      IVWB"2  l/w  1  2  1  2 

2r 

-  Cos6     +  Cos6    ] 


(12A) 
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Hence  expressions  for  B  ,  B~  and  B  „  are 

U   0      y0 


■192   E 
j  =  l 


n  2 


(2j  +  1)! 


*4i  2 

U  J   K.  I 


2(j-3) 


6j   T*2j 


(13A) 


*2    2   3 
where  y   =  u   /ca     and 


E    (2k)!/(k!) 
k=o 


(14A) 


21  _© 
5 


*2 


10 


216  0 
245 


*6 


*3   X15  + 


(15A) 


*2    2,5 
where  0   =0  /ea   and  finally 


5y0 


*?   *2 

3u   0 

"  *2 


*2 


I„  - 


11 


72 

35 


*2  4 

y   0^ 

m*3 


X13  + 


(16A) 


Note  that  equations  (13A) ,  (15A)  and  (16A)  are  expansions  for  which  only  the  first 
few  terms  are  considered.   The  index,  j,  of  equation  (13A)  is  set  equal  to  4. 

The  above  equations  contain  integrals  I   given  by 


exp  [-B4>s  (r  )]  (r  )  X     dr 


(17A) 


where  r   =  r/a. 

The  Pople  expansion  technique  has  the  advantage  that  the  integration  over 
angles  only  has  to  be  carried  out  once.   Hence  the  calculation  of  the  second 
virial  using  a  particular  spherical  potential  reduces  to  the  straightforward 
integrations  of  equation  (17A) ;  the  results  of  which  can  be  tabulated  for  a 
given  spherical  potential. 
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f 
Returning  to  nitrous  oxide:   values  for  the  moments  are  [38] 


0  =   3.0  x  10  26  e.s.u. 


(18A) 


—1  8 
u  =  0. 166  x  10    e. s.u. 


Insertion  of  these  parameters  into  the  expression  for  B  ,  equation  (9A) ,  led 
to  the  graph  shown  as  figure  3b.   Also  included  in  this  figure  are  the  data 
points  of  Hirth  and  Kobe  [30]. 

It  can  be  seen  that  our  prediction  of  the  second  virial  is  not  unsatisfac- 
tory but  the  predicted  values  are  consistently  somewhat  higher.   Possible 

J. 

reasons  for  the  discrepancy  are  1)  neglect  of  the  induction  contribution  to  B  , 
2)  lack  of  convergence  of  the  expansions  since  the  quadrupole  moment  of  nitrous 
oxide  is  relatively  large  and  3)  uncertainty  in  the  m-6-8  parameters.   Unfortunately, 
we  cannot  correct  this  last  factor.   We  have  shown  that  the  m-6-8  parameters 
should  be  selected  from  viscosity  data  taken  over  a  wide  temperature  range  (see 
the  discussion  in  reference  [32])  but  the  data  is  limited  for  nitrous  oxide. 
We  did,  however,  investigate  2)  as  remarked  in  the  next  section. 

Direct  Calculation  of  the  Second  Virial  Coefficient.   There  are  some  drawbacks 
in  the  expansion  (9A) .   For  example,  all  I  integrals  have  to  be  redetermined  if 
the  spherical  potential  is  modified  and,  of  interest  here,  the  arbitrary 
truncation  of  the  expansions  can  introduce  convergence  problems.   Several 
authors  have  investigated  alternative  methods  to  the  expansion  of  the  exponential 
term  in  equation  (8A) .   The  most  direct  is  to  carry  the  full  integration 
numerically  using,  for  example,  the  Monte-carlo  approach.   This  and  other 
integration  schemes  are  discussed  by  Stroud  [39]. 


Steele  and  Sweet  [40]  discuss  a  procedure  in  which  the  potential  function, 
<Kr,  9,  9„  (j))  ,  and  the  radial  distribution' function,  g(r,  6  9„  cf>)  ,  are  expande 
in  spherical  harmonics  and  the  integrations  over  angles  performed  numerically. 


1  e.s.u.  =  3.33564  x  10  U   C  m2 
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The  virial  coefficient  follows  from  a  straightforward  integration  over   r 
intermolecular  separation. 


the 


Our  work  with  nitrous  oxide  here,  and  previously  with  carbon  dioxide  [34], 
suggested  another  integration  technique  which  avoids  the  somewhat  complex  and 
time  consuming  calculations  alluded  to  above,  but  which  can  eliminate  the 
convergence  problems  of  the  expansion  method.   An  outline  of  our  procedure  is  as 
follows. 

Consider  the  potential  of  equations  (3A),  (10A)-(12A)  and  the  coordinate 
system  of  figure  4.   The  two  molecules  are  set  at  fixed   r   apart  at  some  initial 
relative  angular  orientation,  e.g.,  6   =  6   =  (J)  =  0.   The  angles  are  then 
varied  incrementally  (say  in  increments  of  tt/4,  tt/8,  etc.)  over  their  range  of 

tt   c})  =  2tt)  and  the  numerical  value  of  the  angle  dependent 
part  of  the  potential  recorded  for  each  increment.   Due,  however,  to  the 
symmetry  of  the  model,  many  of  the  values  are  identical  so  only  the  different,  e, 
values  and  their  weight,  wt,  (their  fraction  of  the  total)  are  recorded.   The 
second  virial  coefficient  for  a  given  initial  configuration  follows  after 
integrations  over   r   according  to  the  expression 


integration  (0  ,  0 


2ttN 


e 

Z 

i=l 


(sswt) 


1 

1  J  o 


(1 


-4).(r,  e1,  e2,  4>) 

kT 


9 

r  dr 


Z   (sswt) 
i=l 


(19A) 


where 


(sswt) 


(Sin©.,  Sin0n  wt)  . 
1      2     l 


i  =  1,  e 


(20A) 


Details  of  the  calculation  will  be  published  elsewhere.   Values  of   B   for  two 
initial  configurations  at  two  temperatures  were  determined  and  are  shown  as 
triangles  in  figure  3b.   One  can  conclude  that  the  expansion  procedure,  shown 
by  the  curve,  did  indeed  underestimate  slightly  the  virial. 
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It  should  be  stressed  that  the  calculations  of  the  second  virial  coefficient 
are  predictions :   one  could  get  better  agreement  with  experiment,  without 
affecting  significantly  the  viscosity  fit,  by  adjusting  the  potential  parameters. 

Finally,  we  conclude  from  the  calculation  of  the  second  virial  that  our 
assessment  of  about  +  3%  accuracy  in  the  viscosity  data  is  reasonable. 
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Figure  3(a).   Viscosity  coefficient  of  dilute  N„0:   comparison  between 
kinetic  theory  values  and  equation  (31).   3(b).   Plot  of  the  second  virial, 
B,  versus  temperature.   Data,  filled  circles,  from  reference  [30].   The  curve 
is  from  equation  (9A) ,  while  the  triangles  are  results  discussed  in  the  text. 
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analysis,  theoretical  physics  and  chemistry,  logical  de- 
sign and  programming  of  computers  and  computer  sys- 
tems. Short  numerical  tables.  Issued  quarterly.  Annual 
subscription:  Domestic,  $9.00;  Foreign,  $11.25. 

DIMENSIONS/NBS  (formerly  Technical  News  Bulle- 
tin)— This  monthly  magazine  is  published  to  inform 
scientists,  engineers,  businessmen,  industry,  teachers, 
students,  and  consumers  of  the  latest  advances  in 
science  and  technology,  with  primary  emphasis  on  the 
work  at  NBS.  The  magazine  highlights  and  reviews 
such  issues  as  energy  research,  fire  protection,  building 
technology,  metric  conversion,  pollution  abatement, 
health  and  safety,  and  consumer  product  performance. 
In  addition,  it  reports  the  results  of  Bureau  programs 
in  measurement  standards  and  techniques,  properties  of 
matter  and  materials,  engineering  standards  and  serv- 
ices,  instrumentation,   and    automatic    data    processing. 

Annual  subscription :  Domestic,  $12.50;  Foreign,  $15.65. 

NONPERIODICALS 

Monographs — Major  contributions  to  the  technical  liter- 
ature on  various  subjects  related  to  the  Bureau's  scien- 
tific and  technical  activities. 

Handbooks — Recommended  codes  of  engineering  and 
industrial  practice  (including  safety  codes)  developed 
in  cooperation  with  interested  industries,  professional 
organizations,  and  regulatory  bodies. 
Special  Publications — Include  proceedings  of  conferences 
sponsored  by  NBS,  NBS  annual  reports,  and  other 
special  publications  appropriate  to  this  grouping  such 
as  wall  charts,  pocket  cards,  and  bibliographies. 
Applied  Mathematics  Series — Mathematical  tables,  man- 
uals, and  studies  of  special  interest  to  physicists,  engi- 
neers, chemists,  biologists,  mathematicians,  com- 
puter programmers,  and  others  engaged  in  scientific 
and  technical  work. 

National  Standard  Reference  Data  Series — Provides 
quantitative  data  on  the  physical  and  chemical  proper- 
ties of  materials,  compiled  from  the  world's  literature 
and  critically  evaluated.  Developed  under  a  world-wide 
program  coordinated  by  NBS.  Program  under  authority 
of  National  Standard  Data  Act  (Public  Law  90-396). 


NOTE:  At  present  the  principal  publication  outlet  for 
these  data  is  the  Journal  of  Physical  and  Chemical 
Reference  Data  (JPCRD)  published  quarterly  for  NBS 
by  the  American  Chemical  Society  (ACS)  and  the  Amer- 
ican Institute  of  Physics  (AIP).  Subscriptions,  reprints, 
and  supplements  available  from  ACS,  1155  Sixteenth 
St.  N.W.,  Wash.  D.  C.  20056. 

Building  Science  Series — Disseminates  technical  infor- 
mation developed  at  the  Bureau  on  building  materials, 
components,  systems,  and  whole  structures.  The  series 
presents  research  results,  test  methods,  and  perform- 
ance criteria  related  to  the  structural  and  environmental 
functions  and  the  durability  and  safety  characteristics 
of  building  elements  and  systems. 

Technical  Notes — Studies  or  reports  which  are  complete 
in  themselves  but  restrictive  in  their  treatment  of  a 
subject.  Analogous  to  monographs  but  not  so  compre- 
hensive in  scope  or  definitive  in  treatment  of  the  sub- 
ject area.  Often  serve  as  a  vehicle  for  final  reports  of 
work  performed  at  NBS  under  the  sponsorship  of  other 
government  agencies. 

Voluntary  Product  Standards — Developed  under  proce- 
dures published  by  the  Department  of  Commerce  in  Part 
10,  Title  15,  of  the  Code  of  Federal  Regulations.  The 
purpose  of  the  standards  is  to  establish  nationally  rec- 
ognized requirements  for  products,  and  to  provide  all 
concerned  interests  with  a  basis  for  common  under- 
standing of  the  characteristics  of  the  products.  NBS 
administers  this  program  as  a  supplement  to  the  activi- 
ties of  the  private  sector  standardizing  organizations. 
Consumer  Information  Series — Practical  information, 
based  on  NBS  research  and  experience,  covering  areas 
of  interest  to  the  consumer.  Easily  understandable  lang- 
uage and  illustrations  provide  useful  background  knowl- 
edge for  shopping  in  today's  technological  marketplace. 

Order  above  NBS  publications  from:  Superintendent 
of  Documents,  Government  Printing  Office,  Washington, 
D.C.  20402. 

Order  following  NBS  publications— NBSIR's  and  FIPS 
from  the  National  Technical  Information  Services, 
Springfield,  Va.  22161. 

Federal  Information  Processing  Standards  Publications 
(FIPS  PUBS) — Publications  in  this  series  collectively 
constitute  the  Federal  Information  Processing  Stand- 
ards Register.  Register  serves  as  the  official  source  of 
information  in  the  Federal  Government  regarding  stand- 
ards issued  by  NBS  pursuant  to  the  Federal  Property 
and  Administrative  Services  Act  of  1949  as  amended, 
Public  Law  89-306  (79  Stat.  1127),  and  as  implemented 
by  Executive  Order  11717  (38  FR  12315,  dated  May  11, 
1973)  and  Part  6  of  Title  15  CFR  (Code  of  Federal 
Regulations). 

NBS  Interagency  Reports  (NBSIR) — A  special  series  of 
interim  or  final  reports  on  work  performed  by  NBS  for 
outside  sponsors  (both  government  and  non-govern- 
ment). In  general,  initial  distribution  is  handled  by  the 
sponsor;  public  distribution  is  by  the  National  Techni- 
cal Information  Services  (Springfield,  Va.  22161)  in 
paper  copy  or  microfiche  form. 


BIBLIOGRAPHIC  SUBSCRIPTION  SERVICES 


The  following  current-awareness  and  literature-survey 
bibliographies  are  issued  periodically  by  the  Bureau: 
Cryogenic   Data  Center  Current  Awareness  Service.  A 

literature    survey   issued   biweekly.   Annual   subscrip- 
tion: Domestic,  $20.00;  Foreign,  $25.00. 
Liquified  Natural  Gas.  A  literature  survey  issued  quar- 
terly.  Annual   subscription:    $20.00. 


Superconducting  Devices  and  Materials.  A  literature 
survey  issued  quarterly.  Annual  subscription:  $20.00. 
Send  subscription  orders  and  remittances  for  the  pre- 
ceding bibliographic  services  to  National  Bureau  of 
Standards,  Cryogenic  Data  Center  (275.02)  Boulder, 
Colorado  80302. 
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